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Abstract 

Flat FRW perfect fluid cosmologies can be reproduced as particular 
solutions of suitable field theoretical models. Here we investigate the sta- 
bility of perfect fluid model trajectories with respect to sets of trajectories 
of the corresponding field models having generic initial conditions. It is 
shown that the trajectories of barotropic perfect fluid models and those 
of the Chaplygin gas model are stable. The total probability to reach the 
Chaplygin gas regime early enough to achieve a matter dominated stage 
having a realistic duration is calculated for a scalar field model assuming 
a natural measure in the space of initial conditions taken near a cosmo- 
logical singularity. An example is presented of a two-fluid cosmological 
model where the stability is partially absent. 

1 Introduction 

The discovery of the phenomenon of cosmic acceleration 012] has stimulated an 
intense activity directed to the construction of new cosmological models based 
on different microscopical models of dark energy producing this acceleration, in 
particular, those unifying dark matter and dark energy (see, e.g., the reviews 

0CQED 

It is well-known that for isotropic flat Friedmann-Robertson- Walker (FRW) 
cosmological models, for a given dependence of the cosmological scale factor a(t) 
on time t, it is always possible to construct a potential U(<p) for a minimally 
coupled scalar field <p which reproduces this cosmological evolution exactly (see, 
e.g., HHED- A similar statement holds also for tachyon matter (for a detailed 
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analysis and examples see [S]). If one has an exact solution arising in a FRW 
universe filled with a perfect fluid, one can always implicitly (and often explicitly 
as well) find a scalar field or tachyon model which reproduces this evolution 
provided suitable initial conditions are imposed However, it is necessary 

to emphasize that the correspondence between perfect fluid, scalar field and 
tachyon cosmological models has a very limited character and amounts precisely 
to the existence of identical solutions obtained thanks to a special choice of initial 
conditions. If one moves away from these conditions, the dynamics of a scalar 
field FRW model can be very different from that of a corresponding perfect fluid 
model or of a model with tachyon matter. 

In spite of the existence of a vast variety of scalar field, tachyon |5J ^| ^1 
IT^l IT3l mi fc-essence QS] an d other field-theoretical cosmological models, 
models based on perfect fluids still retain their own attractive features due to 
their simplicity and a smaller number of parameters. 

Consider, for example, the Chaplygin gas cosmological model introduced in 
[TC] and further developed in |17l 1181 IT§1 I2U) . This model is the simplest one 
unifying dark matter and dark energy in a FRW background. It is based on the 
equation of state 



where A is a positive constant. In what follows, we shall consider only the 
Chaplygin gas cosmological model with p > \f~A which is equivalent to the 
energy dominance requirement that \p\ < p. The Chaplygin gas with p < \f~A 
was considered elsewhere E3- (While this latter case is obviously not apt to 
the unification of dark energy and dark matter, it could be useful in the context 
of phantom cosmology |2~2"]). 

It was shown in ^H] that the Chaplygin gas cosmological dynamics can be 
reproduced using the scalar field model with the standard kinetic term ^(p^c/)'^ 
and the potential 



(throughout the paper, we use the Landau-Lifshitz sign conventions, in particu- 
lar, the space-time signature (H ), and units in which c = h = 8ttG/3 = 1, 

then A — p 2 min <C 1). This representation is unique up to a constant shift of <f>. 
However, the opposite is not true: not any (and in fact only one) FRW solution 
for a minimally coupled scalar field with the standard kinetic term and the po- 
tential J5J satisfies the relation Thus, there arises the problem of stability 
of dynamical trajectories of the FRW Chaplygin gas model with respect to the 
family of trajectories of the scalar field FRW model with the potential (J2J. 

It has also been noticed that the cosmological model based on Sen's effective 
Lagrangian of tachyon matter [Hj 





L = -V(T) y /l-T ll T<i> 



(3) 
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with the potential V(T) = y/A exactly coincides with the Chaplygin gas model 
This shows that the representation J5J is not the only possible microscop- 
ical model for the Chaplygin gas if the assumption of the standard (minimal) 
kinetic term is dropped. Moreover, remarkably, the tachyon representation of 
the Chaplygin gas has no extra solutions with a time- like gradient: all solutions 
of the model J2J with T^T^ > satisfy the equation of state |QJ. 

It is known that the Chaplygin gas model has serious problems with linear 
perturbations on the FRW background: analysis shows that it has a wrong 
transfer function for matter perturbations and a too low value of high 
multipole temperature anisotropy |241 1251 I26| . Account of non- linear effects, 
while leading to some new features, is not sufficient to save the simple Chaplygin 
gas model [21]. Nevertheless, we believe that in spite of this failure, meaning 
only that this model is too simplified to describe the real Universe, the very idea 
of dark matter and dark energy unification remains alive. Models with a more 
complicated effective equation of state resulting in a lower value for the sound 
velocity at recent times may avoid this difficulty. In particular, such models 
may be based on Lagrangians of the type © but with different dependence on 
TpT' 1 * 28 . Also, non-adiabatic perturbations in the Chaplygin gas |2S] and 
more general models open a new possibility to cure this problem. Thus, it 
seems to us reasonable to share the hope, expressed by Bilic et al (2Zj, that "the 
Chaplygin gas model is not so much wrong as incomplete". For this reason, in 
this paper we restrict ourselves to the FRW background in order to investigate 
if the Chaplygin gas model has problems already at this level, as was recently 
claimed in [8*T] . 

It is tempting to conceive that the appearance of the Chaplygin gas in cos- 
mology has a microscopic origin and represents some remnant of a more funda- 
mental theory (see, for example, while it is difficult to imagine that the 
scalar field model with potential has a special meaning. However, just this 
scalar field representation of the Chaplygin gas model was used in Here a 
numerical study was performed of the problem of stability of the Chaplygin gas 
model trajectories in the framework of a wider set of trajectories of the scalar 
field model This numerical analysis was used to draw the conclusion that 
these trajectories are unstable, which was interpreted as a new and very serious 
problem for the Chaplygin gas cosmology arising already at the level of back- 
ground solutions. As we have already noted above, the scalar field model |2j) is 
neither unique, nor a fundamental representation of matter with the Chaplygin 
gas model equation of state. So, the problem of stability of the Chaplygin gas 
model trajectories with respect to the set of trajectories of this model cannot 
serve as a touchstone for the validity of the Chaplygin gas model as a candidate 
for a model unifying dark energy and dark matter. In particular, there is no 
such problem at all for the tachyon representation © . 

However, the problem of stability of trajectories of perfect fluid models with 
respect to a set of trajectories of underlying field-theoretical models is of interest 
by its own. That is why we study this problem analytically in this paper. It will 
be shown that the Chaplygin gas model cosmological evolution is stable with 
respect to small deviations, considered in the linear approximation. Then we 
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consider other perfect fluid models and their scalar field and tachyon "relatives" 
and show that the linear stability of perfect fluid trajectories is a non-trivial 
phenomenon which does not always occur. 

The method to study stability of perfect fluid trajectories is the following 
one. We introduce a parameter measuring the deviation of a scalar field or 
tachyon model trajectory from the corresponding perfect fluid trajectory and 
then study the time evolution of this parameter in the neighbourhood of the 
perfect fluid trajectory in the linear approximation (namely, to first order in 
the parameter). In Sec. II this method will be applied to different perfect 
fluid models in the regime in which deviations between scalar field and perfect 
fluid trajectories are small and, therefore, can indeed be treated to first order. 
Section III will be devoted to the detailed study of the Chaplygin gas model. 
Here we shall abandon the hypothesis of small linear deviations. The total 
probability to reach the Chaplygin gas regime early enough to achieve a matter 
dominated stage having a realistic duration will be calculated for a scalar field 
model assuming a natural measure in the space of initial conditions taken near 
a cosmological singularity. It will be shown that the Chaplygin gas trajectory 
is rather stable not only in a local but also in a global sense. 



2 Linear analysis of the stability of perfect fluid 
cosmologies 

The first model that we investigate is the Chaplygin gas. We introduce the 
parameter 

x = pp + A, (4) 

which is equal to zero for the Chaplygin gas cosmological trajectory. Now, 
consider the scalar field model with the potential For this model, the 

Einstein equations for a FRW space-time with the zero spatial curvature have 
the form 

*-{%)' ->-$™ < 5 > 

h = -~(p + p), p=^--U(d>). (6) 
Substituting the expressions (J5J and © into the formula one gets 

x=^-U\cp)+A. (7) 

The time derivative of x is 

x = <^> 3 - 2UU'4>, (8) 
where U' = dU /d<j>. Using the Klein-Gordon and Friedmann equations, one has 



-3\%- + U4>-U'. (9) 
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Substituting © into Eq. iJSJ, we obtain 



i = -30 4 yy + £/-0 3 [/'-2[/[/'0. (10) 

Now, using Eq. Q , we can express the time derivative of the scalar field as 

<P = -V2(U 2 - A + x) 1 / 4 

(A \ 1/4 
= -V2 ( j tanh 2 30 + a; J . (11) 

Substituting this expression together with the explicit form of the potential J5J 
into Eq. ((171)1 . we get the expression for the time derivative x as a function of 
the field and the parameter x. Precisely: 

x = -3(Atanh 2 30 + Ax) 

1/2] 1 / 2 

-VA ( cosh 30 H ) + ( — tanh 2 30 + x 

2 \ cosh 30 / \ 4 

„ f— sinh 3 30 / A , » , \ 3/4 

-3V2A ^— - tanh 2 30 + x) 

cosh 2 30 V 4 / 

3\/2A sinh 3 30 / , OJ 1 \ (A l2oJ \ 1/4 ,.. 

'cosh30H — - — tanh 2 30 + x) . (12) 



2 cosh 2 30 V cosh307 \4 

We write down the leading term of the expansion of the right-hand side of Eq. 
Ijl2|l in the neighborhood of the point x = 0, disregarding terms of order higher 
than one: 

3A 1 / 4 (5sinh 2 30 + 6) 
X = 2(^30)3" X + ^ ^ 

The sign of the coefficient of the linear term in x is negative. Therefore, tra- 
jectories of the Chaplygin gas cosmological model, if considered as a subset of 
trajectories of the scalar field model (J2J, are stable. 

Next we calculate the rate of convergence of scalar field model trajectories 
to the Chaplygin gas solution x = 0. For the Chaplygin gas the value of the 
Hubble parameter is 

h = A 1/4 v /cosh30. (14) 
Thus, the relation between the coefficient 7 from the equation 

x = —jx + o(x) (15) 

and the Hubble parameter is the following: 

7_3^, 1 , (16) 



h 2 V cosh 2 30 
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and hence 

7-5 < ~ < 9. (17) 

Thus, the Chaplygin gas behaviour is quickly achieved once \x\ ~ A. 

For the tachyon matter representation J2J of the Chaplygin gas model (V(T) — 
\f~A = const) , the variable x is identically zero, so there is no problem of stability 
at all, as has been pointed out above. 

We now consider another well known perfect fluid model obeying the barotropic 
equation of state 

P = kp. (18) 

A spatially flat FRW universe filled with this fluid expands according to the 
power law 

a(t) = a Q t^m . (19) 

The cosmological dynamics (|19fl can be reproduced in the scalar field model 
with the standard kinetic term and the potential [351 |S] 



U((f>) = [/ o exp(3Vl + fc0). (20) 

In order to measure the deviation of a trajectory of the model (|18[) from that 
of the model l|2U|l , we introduce a variable x as 

x=^-k. (21) 
P 

Following the method described above for the Chaplygin gas model, we come 
to the equation 

x = |v/2[/(l- (k + x)2)(y/l + k - y/l + k + x) 



= - - y/2U(l-k)x + o(x) (22) 

and see that the perfect fluid trajectory is stable in this case, too. 
For the power-law evolution Ijl9|l 



"=VS (23) 

and 

J- (^) 

Thus, we have a convergence of trajectories of the scalar field model to that of 
the perfect fluid model, whose rapidity depends on the value of k. 

It is known that in the case of negative pressure k < 0, the cosmological 
evolution given by Eq. 119|l can be reproduced in the framework of the tachyon 
model © [El ED d with the potential 
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Using the known expression for the tachyon energy density and pressure 

\/\-T 2 



one can easily see that 

x - 

Using the tachyon field equation 



p= -V(T)Vl-t 2 , (27) 
x= -(l-T 2 )-k. (28) 



f + 3VV(l-f 2 ) 3/i f + (l-f 2 )^y- = 0, (29) 

one can show that 

and, hence, the perfect fluid trajectory l|19[) is stable with respect to the tachyon 
model having the potential As far as the case k > is concerned, it was 

shown in [H] that the dynamics of the perfect fluid model can be reproduced in 
the framework of the Born-lnfeld action of the type 



L = W(T)y/T tll T<»-l, (31) 

with the potential 



W{T) = tttWW (32) 



4\/fc 
9(1 + k)T 2 

A similar analysis shows that the trajectory (|19fl is stable also with respect to 
the set of trajectories of the model l|31|) . i|32|l . 

The value of the Hubble parameter for these simple tachyon models is 

(33) 



3vT+l?T 
and, hence, 

1 = ^S. (34) 
h 2VTTI 

Thus, the rate of convergence depends on the value of the parameter k. Note 
that in contrast to formula (|24[l. the relation ? in (|34|) is not bounded above 
and tends to infinity when k — > — 1. 

We now consider a more complicated family of cosmological models. In paper 
jHJ, a two-fluid cosmological model was studied. One of this fluids is merely a 
cosmological constant p = —p = —A, while the second fluid obeys the equation 
of state p = kp, where — 1 < k < 1. The cosmological evolution of the model is 
described by the formula 

/ _ x 2/[3(l+fc)] 

( . 3VA(l+j)A 
a(t) = a smh . (35) 
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It has been shown in |S] that the cosmological dynamics can be repro- 
duced in the scalar field model with the standard kinetic term and the potential 



[/(</>) = A 1 + 



• sinh" 



3VTTk<f> 



(36) 



or in the tachyon model with the potential 
V(T) = — A 



sin 



V / A(l + fe)T 



1 - (1 + k) cos 2 



(37) 



The parameter defining a deviation of the cosmological dynamics of the field- 
theoretical models (j3HJ, (|3Tjl from the two-fluid behaviour i|3"5|) looks as follows: 



P + A 

.X = ; K. 

P- A 



(38) 



For the scalar field model with the potential (|36|l . the time derivative of the 
parameter x in the neighbourhood of the point x = is 



3VA(1 - k) 
2cosh5vSM 



cosh' 



3 V / TTI0 1 



x + o(x) 



(39) 



One can see that in the case k < 0, the coefficient of x in the right-hand side 
of Eq. (|3"§|l is always negative. Hence, the trajectory l|3"5|l is stable. If the 
parameter k is positive, the trajectory Q35JI is stable for 



> 



3VT+T 



arccosh 



I - k 



(40) 



and is unstable when the condition l|4l7|l is not satisfied. Thus, the trajectory 
(|35|l becomes unstable close to the de Sitter point (j> = 0, = 0, a/a = \/A. The 
condition (f^TTf) could be also rewritten as 



1 - k 
1 + fc' 



(41) 



where J7a denotes, as usual, the ratio of the cosmological constant to the general 
energy density in two-fluid model. 

As regards the tachyon model l|3"7|) . using expressions iffi^ . (|27Jl and 
after cumbersome but straightforward calculations we come to the following 
equation: 



3 


2fc(l - k) + {k + l 


2 cos 2 


§VA(l + fc)T 


sin 2 


|v/A(i + fc)r] 


2 sin 


|v/A(l + fc)T 


X 


(1 


-k)-(l + k) 


30S 2 


|VA(l + fc)T 





r + o(a;). 
(42) 
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If k < 0, the denominator of the coefficient of x at the right-hand side of Eq. 
<|42l) is positive while the numerator is negative at 



cos 



or at 



VA(l + fc)T 



> 



1 

<2- 



/l 2k(l-k) 
4 + (yfc + 1) 2 



/l | 2fc(l-fc) 



(fc + 1) 2 



(43) 



(44) 



Thus, the perfect fluid model trajectory (|35|l will be stable if either one of the 
conditions or l|4^jl is satisfied. If 



k < 



4V2 
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(45) 



then trajectory (|35|l is stable for any value of the tachyon field T. 

In the case k > 0, the numerator of the coefficient of x in the left-hand side 
of Eq. l|42Jl is always positive and the denominator is negative if 



cos 



-VA(l + fe)T 



1 



> 



1 



(46) 



Thus, the perfect fluid trajectory is unstable in the vicinity of the de Sitter point 
cos 1^(1 + k)T = 0,T= 0,d/a = VX. 
For the tachyon two-fluid model we have 

VA 



h = 



sin (^y/A(l + k)T 



(47) 



and 



3 ((2fc(l - fc) + (1 + fc)2cos 2 Uy/A(l + k)T\ sin 2 (| v /A(l + /c)T 
^ 2 ^((1 + k) cos 2 (f^l + ^r) -(1-fc))) 

(48) 

and different regimes of convergence are possible depending on values of fc, A 
and T. 



3 Analysis of a large class of initial conditions 
for the Chaplygin gas 

In this section we go beyond the linear analysis, developed in the preceding 
section, of the small deviations of the trajectories of scalar field and tachyon 
cosmological models from the trajectories of the corresponding perfect fluid 
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models. For the most interesting case of the Chaplygin gas, we consider a broad 
variety of initial conditions for the corresponding scalar field model. However, 
in order to obtain definite results, we need first to specify which models are 
compared to each other. For realistic cosmological models describing not only 
the present Universe, but the early Universe too, one has to include radiation. 
So, a possible perfect fluid model is the two component one consisting of the 
Chaplygin gas and radiation (p 7 = p 7 /3); we neglect baryons here. But with 
which scalar field model should we compare it with? If we want to introduce 
a minimally coupled scalar field <f> with the standard kinetic term and some 
potential U(tfi) which models the Chaplygin gas only (so that the comparison 
model is scalar field + radiation), this potential should be very peculiar: it 
becomes strongly divergent at some finite value of 4>- 

Indeed, though in this case we cannot obtain an exact form of this potential, 
we can study its asymptotic behavior for p, p 1 ^> \[A and match it to the 
potential @ valid for p^> p 1 . In the former regime, the Chaplygin gas may be 
well approximated by dust-like matter. So, let us consider a scalar field which 
mimics only dust in a universe filled with dust and radiation. As usual, the 
scalar field potential is 

1 C 

where C is a constant describing the amount of dust. Then 



U =o(p-p) = ^s> (49) 



2 -p + P=^, (50) 



C 
a 
or, 

where the squared Hubble parameter is defined from the Friedmann equation 
which contains both dust and radiation: 

" 2 ° ° (52) 
Substituting (|52l into Eqs. ||SDJ|, (|ST|> . one obtains 

dd> 1 



da ^Ja(a + a eq ) ' 

where 



(53) 



a eq = q (54) 

is the value of the scale factor when dust and radiation have equal energy den- 
sities. Introducing a new variable 

a = a eq sinh 2 x, (55) 
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one can integrate Eq. (|53[) obtaining 



= 2\ + 4>a = 2arcsinh A /— + (p , (56) 



Inverting the relation (j56|) . one has 



a = a eq sinh — - — . (57) 

Substituting into Eq. I|49|l . it gives the potential which as announced before 
strongly diverges at <j) — (f>Q 

= i^sr < 58 > 

On the other hand, if the Chaplygin gas dominates over radiation (p» p 7 ), 
the potential U(4>) should have the form J2J). Let us match the expressions J5J 
and l|58|l at the matter dominated stage p ^> VA, p~^> p 7 . The energy density of 
the Chaplygin gas has the following dependence on a(t): p = y/ A + C 2 /a e where 
the constant C is the same as in Eq. I|52[) . Let f2 7 <C 1 be the radiation energy 
density at the "present" moment a — do, p = po = y A/0.7, a/a = Hq = y/po, 
i.e., at the moment when \ptot\l Ptot has the same value 0.7 as in the standard 
ACDM model with Q m = 0.3, ft A = 0.7. Then 

% = pi - A w 0.43A , D = n 7 p at, , (59) 
a 

= _JL_ = o pa « 1.830- . (60) 
ao 1 + z eq 7 \Jpl~ A 7 

In the region 1 -C <C <^>o , 

[7(0) = I Vie 3 * = ^ e 3 ^" 3 ^ , (61) 

4 a eq 

where the first expression for U(4>) is the asymptote of J2J) for <j>^> 1, while the 
second one is the asymptote of (I58|l for cj>Q — <p ^> 1. Therefore, 



1 1 ^ 

e 3 ^ = -^=«-^f, o «|lnO 7 |+O.87»l . (62) 



Thus, inclusion of radiation drastically changes the scalar field potential from 
having an exponential growth at infinity to having a pole at finite value of <fi. 

If, on the other hand, we choose to model both the Chaplygin gas and ra- 
diation by one scalar field, then U{4>) is regular for finite 0, but the model 
itself becomes artificial since, in particular, it misses the fact that radiation has 
a non-zero entropy. Also, in that case wc would investigate the stability of a 
fluid model different from the original Chaplygin gas one. Note once more that 
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no such ambiguity arises in the case of the tachyon representation of the 
Chaplygin gas. 

So, to present methods able to deal with the case of large deviations between 
trajectories of fluid and field models, and to get quantitative results about the 
probability measure of favourable trajectories, we restrict ourselves to the case 
of the pure Chaplygin gas model compared to the scalar field model with the 
potential (J2J, both applied up to arbitrary large curvatures. It is less complete 
from the physical point of view, but more well defined instead. 

Then, the second step is to decide which solutions of the scalar field model we 
consider as favourable ones. A natural choice for this toy model is to choose all 
those solutions for which the Chaplygin gas attractor has been already reached 
by the moment of the matter-radiation equality in the real Universe z — z eq rj 
3450 (the value from the WMAP data is used here). In particular, it may be 
reached earlier. Note that the redshift z ~ 1 should correspond to the moment 
when p ~ \[A. 

We consider the case x~^> A which is typical at sufficiently early time near 
the singularity. In this regime, the kinetic energy of the scalar field much exceeds 
its potential energy. Then 

/2 1 

a{t) oc t 1/3 , (f> = ±— In- + const . (63) 
3 t 

Thus, for \<j)\ ^> 1, the potential U(<j)) decreases oc exp(3|</>|) oc in the 

course of expansion, if <f> and <f> have opposite signs (in the opposite case U((j)) 
even increases oc t^), i.e., more slowly than the kinetic energy which is oc t . 
Therefore, U generically becomes of the order of the kinetic energy at some 
characteristic moment of time t = tu- 

If \(f>u\ = |0(^c/)| ^ 1 at this moment, the long dust-like 'tracking' regime 
follows with 

a(t) oc i 2/3 , <p — ±^ In ^ + const = =p In a + const, U = ^ . (64) 

The latter regime just represents the behaviour of the Chaplygin gas model 
for large p. Thus, for generic initial conditions defined at some sufficiently large 
redshift z ^> z eq 3450 (the value from the WMAP data jH3] is used here for the 
redshift z eq of the matter-radiation equality), the only necessary and sufficient 
condition for the scalar field model to reach the Chaplygin gas trajectory at a 
redshift exceeding z eq is 

\0u\>\nz eq ^8 (65) 

(recall that 4> is measured in units of (3/87TG) 1 / 2 ). The quantity (j>o where the 
potential U diverges in Eq. ifS!^) is of order of In z eq , too. Also, notice that in the 
case x <C —A, the potential energy is much higher than the kinetic energy. A 
simple analysis of the Klein-Gordon equation shows that in this case the scalar 
field 4> is in the fast rolling regime and, after a Hubble time interval ~ h , it 
has no alternative but to follow the tracking attractor 164JI arriving again at the 
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Chaplygin-like regime. It is clear that, for natural initial conditions defined near 
a cosmological singularity when p is of the order of unity (i.e. , of the Planck value 
in usual units), this condition is satisfied for an open set of initial conditions 
with non-zero measure, being therefore generic. It is instructive to compare the 
value (|65[) to a similar quantity which arises in the theory of chaotic inflation: an 
initial value of a homogeneous inflaton scalar field with a polynomial potential 
at the beginning of inflation which is required to have more than 60 e-folds 
of inflation. For the simplest case U = m 2 <j> 2 /2, the standard expression for 
the number of e-folds leads to \<t>i n \ > (| • 60) 1 / 2 w 9. Thus, condition (|63|) is 
even less restrictive than the one which is assumed in the chaotic inflationary 
scenario. On the other hand, two of the three specific initial conditions (taken 
at the redshift z ~ 10 5 ) studied in do not satisfy iJBSJ. The third one - 
their case (a) - does satisfy l|63)l and embodies a model where the Chaplygin 
gas behaviour is valid from sufficiently early time. In [2], it was considered as 
a bad one for the Chaplygin gas model stability because it does not satisfy the 
final condition p ~ y/A at z = 0. However, for given initial values of (f> and 0, 
the redshift z may not be taken arbitrarily; it should be counted back in time 
from the present moment when z = 0, as has been pointed out above. With this 
re-interpretation and proper normalization of z, this case becomes favorable for 
the Chaplygin gas model stability. 

For illustrative purposes, we have performed some numerical calculations 
integrating Eqs. JSJ) and 10. We have studied the most interesting case of 
the transition to the Chaplygin gas behavior in the regime ip S> y/A. In these 
conditions, the Chaplygin gas behavior is practically indistinguishable from the 
dust like one. We plot the time dependence of the quantity x/h 4 which is equal 
tow + A/h 4 where w = p/p is the equation of state parameter. Thus, it just 
reduces to w in the most interesting regime p ^> \J~A. This quantity is equal to 
1 in the kinetic dominated regime and to — 1 in the potential dominated regime 
with p \[A~. On the other hand, the condition \x/p 2 | <C 1 is just the condition 
that we are close to the Chaplygin gas behavior for any p. 

The "realistic" value of A is given by the condition p min = \f~A~ = 0.7Hq 
where Hq is the Hubble constant. The value Hq = 70 km s _1 Mpc -1 corresponds 
to Hq = 3.5 • 10~ 61 in our units. Therefore, A w 10~ 242 , while natural initial 
conditions have to be placed somewhere near t ~ 1. In order to avoid long 
numerical runs, we propose the following scheme. Let us rescale time to t — 
10 53 t, the energy density to p = 10~ 106 p and x = 10~ 212 a;. This means that 
time is measured in units of ~ 1.56 x 10 10 s, i.e., roughly 100 times less than 
the matter-radiation energy equality moment. Then the rescaled value of A is 
10 212 x 10~ 242 = 10~ 30 . All tildas are omitted below. 

Now let us take different initial conditions for <fi and <fi on the circle h = a/a = 
1 at the initial moment chosen conventionally at to = 1 , (strictly speaking, this 
time moment should be shifted by some amount less than 1 depending on initial 
conditions to get an initial singularity strictly at t = 0, but it is not important 
for t > 1). 

The initial conditions which we consider are the following ones: 
a)0 o = O,[/ o = l,0o = 11.98; 
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b) o = l,Uo = 0.5, O = 11.51; 

c) o = -l,£/ o = 0.5, O = 11.51; 

d) o = 1.342, U = 0.1, 0o = H.21; 

e) o = -1-342, C/ = 0.1, 0o = 11-21; 

f) 0o = 1.378, U Q = 0.05, 0o = 10.98; 

g ) o = -1.378, U = 0.05, 0o = 10.98. 

For all the above cases we have made a run up to t = 1000 to show that the 
Chaplygin gas behavior (\x\/p 2 becomes much less than 1) has been reached by 
t ~ 100, i.e., just when the matter dominated stage begins in the real Universe. 
The dependence oix/h 4 is plotted in Fig. 1. For all the cases we have calculated 
also the value of the scalar field 4>u at the moment tjj which is approximately 
equal to 100 for the cases a)-f). The corresponding values of 4>u are the following 
ones: 

a) 0[/ = 8.416; 

b) 0(7 = 8.422; 

c) <j>u = 8.389; 

d) 0!7 = 8.425; 

e) 0i7 = 8.287; 

f) 0tf = 8.429. 

Apparently, all these values satisfy the criterion QG5|). The case g) represents 
the borderline where the condition l|(j5D is saturated. In this case the moment 
tjj s» 130 and (j>u — 8.015. Thus, it is not possible to take initial values of Uq 
smaller than 0.05 for 0o < 0. 

Therefore, we see that with these "late" initial conditions, it is enough to 
have a sufficiently large 0o or J7(0o) to get a long dust-like stage described by 
the Chaplygin gas. It will be shown below that for "natural" initial conditions 
defined at the Planck scale, this requirement may be further relaxed. In other 
words, the earlier we take an initial moment, the larger is the "good" region in 
the space of initial conditions. 

In Fig. 2 we have plotted all seven cases on the time range smaller than 
1000, displaying some additional details of the time evolution of the value of 
x/h A . It is clearly seen that, for 0o > (cases b),d) and f)), transition from the 
kinetic-dominated stage x/h 4 = w = 1 to the Chaplygin- like matter dominated 
stage x/h A = w ~ goes through a short transient potential-dominated regime 
x/h 4 = w = — 1. 

We now return to the calculation of a measure of initial conditions near a 
cosmological singularity leading to a sufficiently long matter dominated stage 
described by the scalar field model of the Chaplygin gas. To be concrete, now 
we introduce two different probability measures in the space of initial conditions. 
These measures are to be fixed on the basis of some additional (usually physical) 
assumptions. Depending on different assumptions, the probability of a set of 
initial conditions may appear both large and small for the same dynamical 
problem. For instance, in the case of a massive scalar field in a closed FRW 
universe, the probability to have a long quasi-de Sitter (inflationary) stage is 
very small if a natural measure in the space of initial conditions is defined at the 
moment of the maximal expansion of a non-singular solution with a bounce j3fi| . 
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while it is close to unity if initial conditions are defined at the Planck energy 
density of the scalar held before (or just at) the beginning of an expanding 
inflationary stage |37] (see also |3*51l3*§| ). 

We calculate the total probability to have condition l|65[l satisfied for two 
different assumptions about a measure in the space of initial conditions for cj) 
and 4> on the curve p = 1 which corresponds to the Planck energy scale. First, 
we assume the uniform measure in <p, i.e. 

dp, = d(/)/2(f> max . (66) 

For this purpose, we first calculate the total probability of the opposite, 
undesired case when \<t>u\ < 8. This clearly requires U <C (t> 2 on the curve p = 1. 
During the kinetic dominated regime l|63[l . the undesired interval \4>u\ < 8 is 
shifted either in the positive or in the negative direction (depending on the sign 
of (j> in Eq. (|63|l ) keeping its width constant. The maximal value of tfi on the 
p = 1 curve is <\> max — hi(l/v^4)/3 ~ 93. 

Therefore, the total probability of undesired solutions is 

P _ = ipi ~ 9 % , (67) 

while the total probability of the opposite, desired case with the condition l|65|l 
being satisfied is 

P+ = 1 - P- « 91% . (68) 

Now, following |37| . we take the uniform measure on the curve p = 1 itself. 
Namely, we introduce the angular parametrization 

6 2 

y = cos 2 a , U{4>) = sin 2 a (69) 

and assume the measure 

da 

dp = — . (70) 

Z7T 

Hence, U'd(f> = 2 sin a cos ada. Note that for the massive scalar field considered 
in |37], the second choice reduces to the first one for |a| <§; 1 (\4>\ <IC 
our case, the total probability of undesired solutions reads 

1 f da , , 3 f , , r- 1 



P— = TT~ / -T-d<t> = — / d<t>VU = -y/Um.max ■ (71) 
Z7T J d(f> ZTT J TT 

Here, it is taken into account that only small regions around the two points 
a = 0, 7r (with (7 « f , but \(j>\ ^> 1 already) contribute. Ui n , m ax is the 
maximal initial value of the potential at the moment t = ti n corresponding to 
H 2 = p = 1 (ti n = 1/3 for the kinetic dominated regime (|63|l ) for which the 

1,-3/2 



Z, 



eq 



condition l|65l) is not satisfied, so that tjj > t eq ~ H 

The maximal initial value of U on undesired trajectories is saturated for 
tu = teq and with <j> and <fi having the opposite signs. Then U(4>(t)) ~ t~ 2 at 

t = t eq , and it grows as t increases ~ t^^ 2 . Therefore, 

TT. r*j +- 2 +v / 2.-V / 2 1Q-32 /yox 
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and 

P_ ~ 1(T 16 , P+ = l- 0(1(T 16 ) . (73) 

Thus, under the second assumption, undesired trajectories have exceedingly 
small probability. The difference between (|67|l . I|68|) and (|73ll illustrates the 
crucial dependence of total probability values on assumptions made about a 
measure in the space of initial conditions. However, in both cases, the probabil- 
ity to reach the Chaplygin gas trajectory sufficiently early to have the matter 
dominated stage in the Universe completely described by the Chaplygin gas 
model is overwhelming, in sharp contrast with the result of |31) . 

Thus, not only are the Chaplygin gas trajectories stable with respect to those 
of the scalar field model with potential but also the probability to reach 
them sufficiently early is close to unity for natural initial conditions defined near 
a cosmological singularity. 

In conclusion, we have shown that stability of trajectories of perfect fluid 
cosmological models with respect to sets of trajectories of more complicated 
field-theoretical models is a property which is not always present, and sometimes 
it is present only for some pieces of a trajectory under consideration. However, 
the Chaplygin gas cosmological model provides an example of a model where 
entire trajectories are stable. 
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